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1 Introduction. Algebro- Geometrical Solu- 
tions and the Reality Reductions. The topo- 
logical charge. 

The algebro-geometrical (or finite-gap) solutions play a very important role 
in the modern theory of soliton equations with periodic and quasiperiodic 
boundary conditions. The main property characterizing these solutions is the 
following: the wave function of the auxiliary linear operator is meromorphic 
in spectral parameter on the finite part of an algebraic Riemann surface. 
It allows to apply the powerful tools of classical algebraic geometry, and to 
write explicit representations for solutions and wave functions in terms of 
Riemann G-functions of several variables. Usually the algebro-geometrical 
solutions are dense in the space of the periodic ones, and if the nonlinearity 
tends to zero, they degenerate to finite Fourier sums. Therefore their role in 
the soliton theory is rather similar to the role of finite Fourier series for the 
linear PDE's. 

The algebro-geometrical solutions of solition equations were first intro- 
duced in the paper [T] in 1974 dedicated to the integration of the Korteveg- 
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de Vries (KdV) equation with the periodic boundary conditions. Real non- 
singular periodic algebro-geometrical KdV solutions can be characterised by 
the following property: the number of gaps (forbidden bands or zones) in the 
spectrum of the auxiliary Schrodinger operator acting in the Hilbert Space 
L/2(R) is finite. This characterization is not valid for other systems (for ex- 
ample the spectrum of non-self adjoint auxiliary operator for the self-focusing 
Nonlinear Schrodinger equation (NLS) has no gaps for any regular potential), 
but neverthereless the algebro-geometrical solutions are called finite-gap. 

The complete theory of algebro-geometrical KdV solutions involving the 
spectral theory of periodic (quasiperiodic) finite-gap Schrodinger operators 
on the line and its unification with algebraic geometry and analysis on the 
Riemann surfaces, time dynamics and algebro-geometric hamiltonian aspects 
was constructed in the works [21 El IH OH EH El (see in the survey article jH] ) 
and in the works [HIGH]- Complex finite-gap solutions of the NLS and Sine- 
Gordon (SG) equations were constructed in [H1H2]- Finite-gap solutions of 
the 2 + 1 Kadomtsev-Petviashvili (KP) eqation were consturcted in J3| on a 
basis of a purely algebraic formulation of the finite-gap approach developed in 
this paper. The algebro-geometric spectral theory of the stationary periodic 
2D Schrodinger operator was developed in the works ^2JE3- In these papers 
it was shown, that the proper analog of the finite-gap constraint in the 2- 
dimensional case is the following: the restriction of the Bloch variety to one 
energy level is algebraic, and the spectral data is "collected" at this energy 
level only. Such potentials can be treated as "integrable" at one energy. 
It means, that the wave functions for this energy and arbitrary complex 
quaismomenta can be written explicitly in terms of Riemann B-functions, 
but the information about the spectral properties at other energies is very 
restricted. 

The algebro-geometrical approach allows to obtain explicit compact for- 
mulas for solutions. But the B-fucntions of several variables are very compli- 
cated from the analytic point of view; therefore the existence of such formulas 
does not lead sometimes to the simple solution of the problem of selecting 
physically or geometrically meaningful solutions among generic complex ones. 
For example, it very easy to describe reality conditions for the KdV, defo- 
cusing NLS or Sinh-Gordon equation 

— =4sinhu(£,77). 

(in all these cases the auxiliary linear operators are self-adjoint). However, 
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the reality reductions for the self-focusing NLS or the Sine-Gordon equation 
found in ^B] are rather nontrivial. It was shown in ^5] that analogous 
reductions are responsible for the absence of magnetic field it the theory of 
the 2D Schrodinger operator. A good review of reality reductions can be 
found if H7|. 

The Sine-Gordon equation 



derived in the XIX century geometry in the light-cone representation 



x = 2(£ + 77), t = 2(£-7j), dt = 2d x + 2dt, d v = 2d x -2d t . (3) 



is one the most important systems of the soliton theory. One it's distin- 
guished properties is the folowing. Usually if we fix the spectral curve (or 
equvalently all local conservation laws), then the level set is a connected Ja- 
coby torus. It is true for KdV, defocusing and self-focusing NLS, the both 
real reductions of the KP equation. But in the Sine-Gordon case the level 
set is a union of 2 m real Jacoby tori, where m denotes the number of pairs 
of negative branch points. As a corollary, the result of averaging a generic 
functional depends not only on the spectral curve, but also on the compo- 
nent. This property is also essential if we calculate the action- angle variables 
(this problem was studied in [T%1 ITUj). 

The averaged densities of conservation laws arose naturally as hamilto- 
nians of the Whitham equations. These equations describe the evolution of 
parametrs of the asymptotic solutions in the form of a slowly varying N- 
phase wavetrains. The Whitham equations for the Sine-Gordon were derived 
in the paper [2U] , but the connectivity components of the level sets were not 
discussed. Let us point out, that the Whitham equations are hyperbolic only 
in the "stable" case, when all branch point are negative and the number of 
components is always greater than 1. 

It is possible to extract from the calculations of j20] , that for all classical 
local conservation laws except the topological charge the result of averaging 
is the same for all components, but this property was never pointed out 
explicitly. In Section 0] we give a simple explanation of it. 



Utt ~ u xx + sinw(x,t) = 



(1) 



= 4 sin w, u = u(£,rj), 



(2) 



where 
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Remark. The number of real tori was calculated in ^H], their character- 
ization in terms of the complex Jacoby torus was obtained in |21j . see also 
|22j . An elementary description of these tori was obtained by the authors in 
[231211, see Section H below. 

Another interesting property of the Sine-Gordon equation is the following. 
In the theory of constant negative curvature surfaces u(£, rj) denotes the angle 
between asymptotic lines. In the theory of Josephson jucntions u has the 
meaning of a phase. In both situations it is defined up to an additive constant 
c = 27m. Therefore it is natural to call a solution u(x, t) space-periodic with 
the period T if exp{iu(x + T)} = exp{iu(x)}. The integer n such, that 
u(x + T) = u(x) + 27m is called Topological charge, the real number 
n = n/T denotes the Density of topological charge. These quantities are 
the "most stable" conservation laws surviving all real periodic nonintegrable 
perturbations. The density of topological charge can be naturally extended 
to all nonsingular quasiperiodic solutions using the formula 

n = lim [u(x + T) — u(x)]/2ttT. 

T — >oo 

One of the most natural problems is the following: How to calculate the 
topological charge of real periodic SG solutions in terms of the 
Spectral Data: the Riemann surface and the divisor? It turned out, 
that nobody succeeded to extract the answer from the explicit description of 
the real components described in terms of the Jacoby torus and B-functional 
formulas, found in the works [2UI22]- An attempt has been made almost 20 
years ago in the work JH] to solve this problem. For this purpose the so-called 
"Algebro- Topological" approach was developed in [TH], and an interesting 
formula has been proposed. However, as it was pointed out in [TH], the idea 
of the proof presented in ^H] works only if the spectral curve is sufficiently 
close to a degenerate one. A complete solution using a new development of 
the main idea of the work [IS] was obtained only recently by the authors in 
the papers [22] (for the so-called stable curves) and [21] for generic ones. Our 
proof is based on a new effective description of the real components in terms 
of divisors. We present a summary of these results in Section [21 

Remark. Formally we can define topological charge for any non-singular 
complex solution. But if the solution has poles, this quantity became ill- 
defined, and no natural non-singularity condition in terms of the spectral 
data for complex potentials is known. It explains why only real solutions are 
discussed. 
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2 Finite-gap Sine-Gordon solutions. 



The "soliton-type" solutions of Sine-Gordon were found already in the XIX 
century using substitution discovered by Bianchi and S.Lie. Now this method 
is called "Backlund transformation" . 

In the early 1970's G.Lamb found out (see that the the Sine-Gordon 
equation is analogous to KdV in the following sence: the so-called "inverse 
scattering method" can be applied to it. The modern approach was started 
in the work The critical point of this approach is the following zero- 

curvature representation for SG 



* 3 



■(E/ + V)tf, % = -(U -V)V, 



(4) 



where 



U 



U(\,x,t) 

V : 



i(u x + u t ) 
-A 



V(X,x,t) 



-i(u x + u t ) 

-\e lu 
r iu 0. 



(5) 
(6) 



To construct generic complex finite- gap SG solutions assume, that we 
have the following spectral data: 

1. A nonsingular hyperelliptic Riemann surface T [/i 2 = R(X)], where 
R(\) = (A - E k ), such that E = and E t ^ Ej, £,j = 0,1,..., 2g. 

k=0 

It has exactly 2g + 2 branch points E = 0, E\, . . . , E 2g , oo; the genus 
of T is equal to g. A point 7 6 T is by definition a pair of complex 
numbers 7 = (A,/i) such, that /1 2 = R(X). 

2. A divisor D of degree g, i.e. set (or formal sum) of g points D = 
71 + • • • + lg- 

For generic data T, D there exists an unique two-component "Baker- 
Akhiezer" vector-function ^(7, x, t) such that 

1. For fixed (x,t) the function ^(7, x, t) is meromorphic in the variable 
7 G T outside the points 0, 00 and has at most 1-st order poles at the 
divisor points 7^, k — 1, . . . , g. 
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2. \?(7, x, t) has essential singularities at the points 0, oo with the follow- 
ing asymptotics: 



i x — t 



*<^ ( > = ( iAt&HOM ) ^ as A — » 0, (8) 

where <j>i(x,t), (f>2(x,t) are some functions of the variables x,t . 

Denote the divisor of zeroes of the first component ip\ by D(x,t) = 
J2jjj(x,t). The vector-function ^(j,x,t) satisfies to the zero-curvature 
equations (@| with potential given by the formula 

„(*,*)=< In ( 9 ) 
(pi{x,t) 

and the function u(x,t) solves the SG equation (j2J. In terms of the divisor 
of zeroes D(x,t) the potenital u(x,t) can be easily written explicitly: 

p iii(x,t) _ £^ Mq) 




The x and t dynamics of the divisor D(x,t) can be described it terms of 
Dubrovin equations (see |23|)- From (fTUj) it follows, that the potential u(x, t) 
is singular at the point (xq, to) if and only if one of the points jk(%a, to) coin- 
cides with or oo. Collisions of 2 or more divisor points result in singularities 
in the solutions of Dubrovin equations, but all symmetric combinations of 
the divisor points including the potential u(x,t) remain non-singular. 

Remark. The first component ipi of the Baker- Akhiezer vector-function 
\I/ is a partial case of the general scalar two-point Baker- Akhiezer functions 
satisfying to the second order linear Schrodinger equation L\^\ = invented 
in the work ^3]. Here we have 
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—A = log 4>\. In this special case our surface V is hyperelliptic, the selected 
points coincide with 0, oo and corresponding local parameters are chosen as 
vA, 1/VA. The second function xp 2 satisfies to similar equation L 2 ip 2 — 
where the operator L 2 is obtained from L\ by the so-called Laplace Transfor- 
mations and vice versa. It is known that cyclic Laplace Chains of the length 
2 lead to the SG equation in the general complex case (see the work [2*7]). 

Lemma 1 Assume, that the spectral data (T,D), D = ji + . . . ,j g satisfy 
the following reality constraints: 

1. The set of branch points of the Riemann surface T contains real points 
or complex conjugate pairs of points only, and all real branch points are 
nonpositive. Without loss of generality we may assume, that the first 
2m + 1 of them E Q = 0, E 1 , E 2 ,. . . , E 2m are real and > E\ > E 2 > 
. . . > E 2m , and for k > m we have E 2 ^ = E 2 k-i, Imi?2fe 7^ 0. Then 
the spectral curve F is real and admits oZ 2 x Z 2 - group of involutions 
generated by the standard holomorphic involution 

a:T^T: a(\, fi) = (A, —fi), (11) 

transposing the sheets, and by the antiholomorphic involution 

r:r^r: r(A, fi) = (A, fi), (12) 

(complex conjugation). 

2. There exists a meromophic differential Q with two simple poles at the 
points 0, oo such, that the zeroes ofQ are located at the points D + tD. 
We shall call such divisors admissible. Without loss of generality we 
can normalize Q to be real: tQ = Q. 

Then the potential u(x, t) is real and non-singular. 

Remark The constraints on the spectral curve we found in ^2]- The 
characterisation of admissible divisors and the proof of non-singularity were 
obatined in the work |16j . 

It is easy to check (see, for example [2*3] ) 

Lemma 2 If D is admissible divisor, then the divisor of zeroes D(x,t) = 
7i (x, t) + . . . + 7 g (x, t) is admissible for all x, t 
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Assume, that we have a real meromorphic differential Q with exactly 2 
simple poles at the points 0, oo. We can associate an admissible divisor D 
to it if and only if 2g zeroes of Q can be can be presented as an union of 2 
subsets D and D\ such, that D\ = tD. It is possible if and only if all real 
roots of Q have even multiplicity. 

Definition 1 Assume that the spectral curve T satisfy the reality condini- 
tions formulated above. A real meromorphic differential Q onT with exactly 
2 simple poles at the points 0, oo is called admissible if the multiplicity of 
all real roots is even. 

Without loss of generality we may assume that the residued of Q at the 
points and oo are equal to +1 and —1 respectively. Then we have 



where P 9 _i(A) is a polynomial of degree at most g — 1 with real coefficients. 
Let us call the polynomial P 9 _i(A) admissible if Q is an admissible differ- 
ential. 

Denote the map associating an admissible polynomial with an admissible 
divisor by II : D —>■ P g _i(X). The inverse map II -1 is multivalued. If P g _i(A) 
is an admissible polynomial, then the equation 



has 2g + 1 roots. One of them is the point (0, 0), other 2g roots form g pairs. 
To define an admissible divisor we have to choose one point from each pair: 
therefore we have at most 2 9 possibilities depending on the number of real 
points. In the generic case we have no real roots, so the number is equal to 
2 9 in this case. 

The description of connected components suggested by the authors in 
|2"3*1 121] is based on the following simple observation: 
Consider a pair of functions of the real variable A: 




(13) 




(14) 



f±W = ± 



(15) 



A 
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We assume, that the functions /±(A) are defined only for such A ^ that 
P(A) > 0, i.e. these functions are real-valued. Let us draw the graphs of 
these functions y = /+(A), y = /-(A), and fill in the following domains: 

\>oy>^. (16) 




Figl. 



Lemma 3 The polynomial P g _i(A) is admissible if and only if the graph 
of P g _i(A) does not cross the black open domains (but it can touch their 
boundaries). 

Now we can associate with each admissible polynomial a Topological 
Type, i.e. a collection of m numbers s k = ±1, k — 1, . . . , m defined by the 
following rule: 

s k = 1 if P g -i(\) > /+(A) as E 2k < A < E 2k -i, 

s k = -1 if P g _x(A) < /_(A) as E 2k < A < En-i. 

Using simple algebraic estimates (see |23j ) it is easy to check: 

Lemma 4 Pac/i topological type is presented by the convex subset in the space 
of polynomials P g -\- Each of these 2 m connected components is non-empty 
and depends continuously on the branch points E\,. . . ,E 2g . 
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The x and t dynamics of admissible divisors may be rather non-trivial (a 
number of numeric experiments are discussed in |22j). In contrast with the 
self-ajoint case (KdV, defocusing NLS, sinh-Gordon) the admissible position 
of a given divisor point depends on the position of other ones (the admissi- 
bility constraints are non-local). The trajectory of a single divisor point may 
be not well-defined because divisor points may collide and a small variation 
of parameters near the collision point results in a bifurcation of trajectories. 
If the potential u(x,t) is periodic with the period T, the corresponding di- 
visor is generically periodic only modulo permutations. But the following 
properties can be easily proved using the characterization formulated above. 

Lemma 5 Projections of the points of admissible divisors to the X-plane 
could not lie in the segments [Ei,E ], [E 3 ,E 2 ],. . . ,(—oc, E 2m \. Moreover, it 
is possible to choose an open neighborhood of these segments such, that the 
projections of these points could not lie in it. 

Lemma 6 Assume, that the projection of a point 7 S = (A s , fi s ) of an admis- 
sible divisor to the X-plane lie in the segments [E 2 k, E 2 k-i}; then fi s < if 
Sk = 1 and n s > if Sk = — 1. 



3 The topological charge of real solutions 



It is well-known (see, for example j2HI22]) that after the Abel transform the 
x-dynamics of Sine-Gordon corresponds to the motion along the straight line 



X = X + Ux, 



(17) 



in the Jacoby torus, where U denotes the noramalized vector of 6-periods of 
the quasimomentum differential dp. 

1 



= -- 



2tt 



dp. 



(18) 



Let us recall, that by definition dp is a meromorphic differential with zero 
a-periods and two second order poles at 0, oo such, that 



dp 



8V\ 
1 

8\V\ 



A 



A 



dX as A — > oo, 



dX as A — » 0. 



(19) 
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The real part of the Jacoby torus is isomorphic to the factor IR n /Z n , therefore 
we can apply the following simple analytic lemma: 

Lemma 7 Letu(X), X e IR n be a smooth function inW 1 such, that exp(iu(X)) 
is single-valued on the torus M n /Z n ; i.e. exp(iu(X + N)) = exp(iu(X)) for 
any integer vector N. Denote by u(x) restriction ofu(X) to the straight line 
X = Xq + x ■ U. Then the density of topological charge n = lim [u(x + T) — 

T — >oo 

u(x)}/2nT is well-defined; it does not depend on the point Xq and can be 
expressed by the following formula: 

n 

n = Y j n k U\ (20) 

k=l 

where U = (U 1 , U 2 , . . . , U n ), and are topological charges along the basic 
cycles Ah, k = 1, . . . , n: 

u(X\ X 2 ,...,X k + l,...,X n )- u(X\ X 2 , . . . , X\ . . . , X n ) = 2vm fc . 

(21) 

The main step is the calculation of the the charges along the basic cycles 
Ak- To obtain a convenient answer the choice of canonical basis is very 
important. For stable surfaces (surfaces without complex branch points) 
the proper system of a-cycles was constructed in [18J. 




Fig 2. 
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We work with the following system of cuts [E 1 , E ], [E 3 , E 2 ],. ■ ■ , (—00, E 2g \. 
Let us denote the sheet containing the line A 6 1, A > /i > by G + and 
the second sheet by G_ (see Fig 2). 

For the generic real curves with 2m negative branch points E±, ... , E 2m 
and g — m complex adjoint pairs E 2 j-i, E 2 j = E 2 j-i, j — m + 1, . . . , g we 
choose first m a-cycles exactly as in the stable case. The last g — m a-cycles 
should be chosen as coverings on V over the path on the A-plane connecting 
the points -E^j-i and E 2 j for j > m. All these path should not meet each 
other. All of them should cross positive part of the line A > in one point Kj 
(i.e. in two points of Riemann surface) such that K m+ \ < K m+2 < . . . < K g . 
This basis depends on the order of the complex conjugate pairs only (see 
Fig 3). 




Fig 3. 

We assume also that the cycles b\, b 2 , . . . , b m are the ovals lying over the 
segments [E 2 ,E^\, [£4,^3], [#2m, #2m-i] respectively. 

For the complex conjugate pairs of branch points E 2 j, E 2 j-i, j > m, we 
make the cuts connecting each of these points with ioo or — ioo such that the 
cuts do not intersect the real line. As before we denote the sheet containing 
the line A 6 R, A>0/i>0by G + and the second sheet by G_. The orien- 
tation of the cycles bi, I — 1, . . . ,m coincides with the standard orientation 
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of the negative semi-line R_ at the sheet G + ; they are opposite at the sheet 
G_. 

Consider a basic cycle Ak on the real component of Jacoby torus, repre- 
sented by the closed curve B k on this torus. The image of this cycle in the 
surface V (i.e. image of it under the inverse Abel map) is a closed oriented 
curve C k (may be consisting from several connected components) homologi- 
cal to the cycle a k G Hi(T, Z). From LemmaElit follows, that this curve does 
not touch the closed segments on the real line [—00, £?2 m ], . . . , [E 3 , E 2 ], [Ex, 0], 
and if it crosses the negative semi-line at the inteval [E21, E21-1], then it hap- 
pens in the sheet G + if (— l) l ~ 1 si > and in the sheet G- if (— l)' _1 s; < 0. 
From Formula it follows, that the topological charge along the cycle A k 
is given by 

n k = C k oR_, (22) 

where C k is the projection of C k to the A-plane, and o denotes the intersection 
index of curves on the A-plane. Taking into account that any intercetion with 
the negative semi-line is the intersection with one of the 6-cycles we obtain: 

Lemma 8 The topological charge along the cycle A k on the real torus in the 
Jacobian variety is equal to the intersection index of the corresponding curve 
C k with the Counting Cycle of the real component on the Riemann 
surface T 

n k = C k o (f^i-iy-'sibA . (23) 

Let us mention that the cycle C k is homotopic to a k on V 
C k o bi = 5 k i and 

J n k — ( — l) fe for k < m 
1 n k = for k > m. 

Combining (|2U|) with (|24|) . we obtain our main result: 

Theorem 1 For any real solution the density of topological charge along the 
variable x is given by the formula 

1 m 

n = ~ 2 ^(-l)*- 1 ^, (25) 

k=l 

where the components U k of the vector U are defined by Mfy) . 



Therefore 
(24) 
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4 Averaging of the local conservation laws 
and topological types 

For a generic functional 

F l u ] = J fi e%U i e ~ %U -. u x, Uxx, -..,ut, u tx , utxx, ■ ■ -)dx (26) 

the result of averaging over a real componet of the Jacoby torus depends on 
the topolgical type. For example by avearging of u x we get the topological 
charge. But if we average the Sine-Gordon Hamiltonian 



H[u) 



uf u 2 



2 + y + (1 - cos(u)) 



dx (27) 



the result is the same for all components. Therefore it is natural to formulate 
the following: 

Problem. How to characterize the functional such, that the averaging 
does not depend on the topological type? 

A complete classification seems to be a complicated problem. But it is 
rather easy to formulate a natural sufficient condition. 

If f(e lu , e~ lu , u x , u xx , . . . , u t , Utx, Utxx, • • •) is a polynomial of its arguments 
then it can be written as a symmetric rational function of the divisor coor- 
dinates / = f(Xx(x, t),..., X g (x, t),fii(x, t), . . . , fig (X, 

Lemma 9 Consider the following form on the T x T x . . . T 

g times 

W f = /(A a , ...A (/ ./, ; fi g ) ]J(X k - Ai) f[ — (28) 

k<l k=l ^ k 

If we fix all divisor points except 71 we obtain a meromorphic differntial on 
T with poles only at and 00. Assume that the residues of this differential 
are equal to identically in 72,. . .j g . Then for a given spectral curve the 
averaging does not depend on the topological type 

Proof. Following (2H| we can calculate the averaging of F[u] using the for- 
mula: 

F\u\ =C j> ... j>W f (29) 
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where C is the normalization constant, and all cycles a*, have the proper 
topological type. The last condition was not discussed in [20 . The defor- 
mations of the integration cycles do not affect the integral if we do not cross 
0. If the condition of Lemma El is fulfilled, then the we can move the inte- 
gration path through and "reach" any topological type without changing 
the integral. 

As a corollary we immediately get the following statement: 

Theorem 2 Assume that the functional F[u] is invariant with respect to the 
following symmetries: a k : (X k , fa) — > (X k ,-fa), cr k : (Xi, fa) -> (X h fa) for 
I ^ k. Then the result of averaging does not depend on the topological type. 

Important example. Consider the densities of "higher SG hamiltoni- 
ans" defined as odd expansion coefficients of the function 



^x{X,X,t) 

ty(X,x,t) 

at the points and oo. Dircet calculation shows, that 



(30) 



1 



where 



* 4 



* (X-X 1 (x,t))...(X-X g (x,t)Y 

% L= . \p + XQ r >(\)] [-Ai(x, t), . . . (-X g (x, t))} 
* *A(A - X 1 (x,t)) ... (A - X g (x,t)) y /E 1 . . .E 2g ' 



(31) 

(32) 
(33) 



Q^(X), Q V (X) are polynomials of degree g — 1 determined by the following 
conditions: 

Qt(X k (x,t))=fa(x,t), Q^X k (x,t)) = ^p^-. (34) 

Xk[x,t) 

We see that the odd part of ^x^ 1 depends on Xk(x, t) but not fa(x,t). 
Therefore the averaging of the "higher SG hamiltonians" does not depend 
on the topological type. 
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